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Abstract 

In this paper, we present a unified physical picture for the electrostatic at- 
traction between two coupled planar Wigner crystals at finite (but below their 
melting) temperature. At very low temperatures, we find a new regime where 
the attraction, arising from the long-wavelength excitation of the plasmon 
mode, scales with the interplanar distance d as dT 2 . At higher tempera- 
tures, our calculation agrees with known results. Furthermore, we analyze the 
temperature dependence of the short-ranged attraction arising from "struc- 
tural" correlations and argue that thermal fluctuations drastically reduce the 

strength of this attraction. 
61.20.Qg, 68.65.+g, 05.70.Np 
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I. INTRODUCTION 



Electrostatic interactions play an important role in a system of charged macroions in 
an aqueous solution of neutralizing counterions 0. The macroions may be charged mem- 
branes, stiff polyelectrolytes such as DNA, or charged colloidal particles. Recently, there 
has been a great interest in understanding the attraction arising from correlations between 
highly-charged macroions as evidenced in experiments |2j and in simulations M. This at- 
traction cannot be explained by the standard Poisson-Boltzmann (PB) treatment, even for 
an idealized system of two highly charged planar surfaces, with counterions distributed be- 
tween them, since PB, being a mean- field theory, neglects correlations. Indeed, it has been 
proven recently that PB theory predicts only repulsions between two likely-charged objects 
0. Recall that the PB solution |I| for a single charged surface with charge density en 
- where e is the elementary charge and n the areal density - immersed in a solution of 
neutralizing counterions of valence Z, predicts a length scale Xqc — 1/ (2iil B Zn) (where 
Ib = ~ 7 A is the Bjerrum length below which electrostatics dominates the thermal 
energy in an aqueous solution of dielectric constant e = 80 (H2O), ks is the Boltzmann 
constant, and T is the temperature.) Physically, this Gouy-Chapman length Xqc defines a 
sheath near the charged surface within which most of the counterions are confined |J . For a 
moderately charged surface of n ~ 1/100 A~ 2 , Xqc is of the order of few angstroms, and for 
highly charged surfaces and multivalent counterions Z > 1, we have Xqc < Ib, signaling the 
breakdown of PB theory. In this limit, fluctuations and correlations about the mean- field 
potential become so large that the solution to the PB equation is no longer valid ||. 

To account for the attraction arising from correlations, two distinct approaches have 
been proposed. The first approach, based on charge fluctuations, treats the "condensed" 
counterion fluctuations in the Gaussian approximation. This theory predicts a long-ranged 
attraction which vanishes as T — > 0. In the other approach based on "structural" 
correlations first proposed by Rouzina and Bloomfield ||, the attraction comes from the 
ground state configuration of the "condensed" counterions. Indeed, at low temperature, the 



2 



"condensed" counterions crystallize on the charged surface to form a 2D Wigner crystal. 
When brought together, the counterions of two Wigner crystals correlate themselves to 
minimize the electrostatic energy. These staggered Wigner crystals attract each other via a 
short-ranged force that is strongest at T = 0. Although the physical origin of the attraction 
is clear in each approach, the relationship between them remains somewhat obscure and 
their results in the T — > limit are somewhat contradictory. Therefore, it is desirable to 
formulate a unified approach which captures the physics of both mechanisms. 

To this end, we attempt in this paper to develop a detailed physical picture of the 
electrostatic interaction at finite temperature between two planar Wigner crystals in the 
strong Coulomb coupling limit. Since correlation effects are essentially two-dimensional, 
we consider a model system composed of two uniformly charged planes a distance d apart, 
each having a charge density en. Confined on the surfaces are negative point-like mobile 
charges of magnitude e. In order to understand correlation effects that are not captured 
by PB theory, we assume that the charges form a system of interacting Wigner crystals 
(see Fig. [I]). In particular, we compute the electrostatic attraction between the two layers 
by explicitly taking into account both correlated fluctuations and "structural" correlations. 
(By "structural" correlations, we mean the residual ground state spatial correlations which 
remain at finite temperature.) In the former case, we obtain a long-ranged force (~ l/<i 3 ), 
which agrees with the result based on Debye-Hiickel (Gaussian) approximation 0. For the 
latter, a simple expression for the short-ranged force is derived, which shows that thermal 
fluctuations reduce its range, and which in the T — > limit agrees with the known expo- 
nentially decaying result [§,[|. Furthermore, we argue that at zero temperature, there must 
also be a long-ranged force derived from the quantum fluctuations of the plasmons ||, in 
addition to this zero-temperature exponentially decaying force. At very low temperatures 
in the quantum regime, we obtain a new length scale \l (to be defined below), where the 
attraction scales like d~ 2 when d < Xl. 

It should be pointed out that we have assumed a uniform charge distribution on the 
surface of the charged plates in our model for electrostatic attraction, mediated by the 
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"condensed" counterions. This assumption of a uniform neutralizing background may not 
be a good approximation to real experimental settings, since charges on macroion surfaces 
are discrete. For monovalent counterions, they tend to bind to the charges on the surface and 
form dipolar molecules. Therefore, the ground state for this system may not be a Wigner 
crystal, which relies on mutual repulsion among charges for its stability, and short-ranged 
effects are likely to be important. However, for polyvalent counterions, a Wigner crystal is 
likely to form since each counterion does not bind to a particular charge on the surface, and 
a uniform background may be more appropriate. The detailed structure of the ground state 
as determined by short-ranged effects and valences will be the subject for another study. 

Another point worth mentioning concerns the ordering of 2D solids which exhibit quasi- 
long-range-order (QLRO) fL0| . It is well-known that a true long-range order is impossible 
for 2D systems with continuous symmetries. For a 2D solid, which may be described by 
continuum elasticity theory with nonzero long-wavelength elastic constants, the Fourier com- 
ponents of the density function n(r) = J2g ^G( r )e* G ' r average (thermally) out to zero for 
a nonzero reciprocal lattice vector G, i.e. (%(r)) = (e* G ' u ^ r ^) = 0, where u(r) are the 
displacements of the particles from their equilibrium positions, while the correlation function 
decays algebraically to zero: (nc(r)fiQ(O)) ~ t -vg(X) w jth r]c(T) = kB ^^ l +x)" > > wnere A 4 
and A are Lame elastic constants. This slow power-law decay of the correlation function 
is very different from the exponential decay one would expect in a liquid. Hence the term 
QLRO. For a single 2D Wigner crystal, QLRO implies that the thermal average of the 
electrostatic potential at a distance d above the plane is zero at any non-zero temperature, 
in contrast to a perfectly ordered lattice (T = 0) where the electrostatic potential decays 
exponentially with d. This may lead to the conclusion that at finite temperatures the short- 
ranged force between two coupled Wigner crystals should likewise be zero. As we show 
below, this is not the case because the susceptibility, which measures the linear response of 
a 2D lattice to an external potential, nevertheless diverges at the reciprocal lattice vectors 
as in 3D solids HIT . 



This paper is organized as follows. In Sec. [Dj we derive an effective Hamiltonian which 



describes two interacting planar Wigner crystals starting from the zero temperature ground 
state. The total pressure is then decomposed into a long-ranged and a short-ranged com- 
ponent, which are evaluated in Sec. [11 A| and [11 B| , respectively, and a detailed discussion of 
our results is presented in Sec. [11 C . In Sec. |T|, we present an argument for a long-ranged 



attractive force arising from the zero-point fluctuations at zero temperature. In addition, 
we use the Bose-Einstein distribution to calculate the attractive long-ranged pressure in the 
quantum regimes. 

II. EFFECTIVE HAMILTONIAN AND PRESSURE 

We start with the Hamiltonian for two interacting Wigner crystals: TC = Ho + Hint- 
Here, 7io is the elastic Hamiltonian for two isolated Wigner crystals |12| 



W> = \ E / (0 n ^(q) ^(qK^-q), (i) 

where (3~ l = ksT, u^(q) is the Fourier transform of the in-plane displacement field of the 
charges in the iih layer (i = A or B), H a p(q) = 27r ^" 2 P^p + ^P^p Q 2 1S the dynamical 
matrix, ji « 0.245 n 3 ^ 2 l B is the shear modulus |L3| in units of k B T, P^g = q a Qp/Q 2 and 



Pad = $af3 — QaQp/Q 2 are longitudinal and transverse projection operator, respectively. Here, 
Greek indices indicate Cartesian components. Hi nt is the electrostatic interaction between 
the two layers: 

= Is /cftcdV^W-"]^-"!, (2) 
J ^(x-x') 2 + rf 2 

where Pi(x) is the number density of charges in the iih layer. In order to capture the long- 
wavelength coupling as well as discrete lattice effects which are essential for our discussions 
on the short-ranged force, we employ a method, similar to that in Ref. [H], which allows 



us to derive an effective Hamiltonian that is valid in the elastic regime where the density 
fluctuations are slowly varying in space, i.e. V • u' ! '(x) 1, but |u A (x) — u B (x)| need not 
be small compared to the lattice constant a. 
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Let us introduce a slowly varying field for each layer: 

4®(x)=z a -u®$*(x)], (3) 

where the displacement field u^(x) is defined in such a way that it has no Fourier compo- 
nents outside of the Brillouin Zone (BZ). Then, the density p«(x) can be written as: 

Pi(x) = £<5 2 (R - ?°(*)) det[a a ^(x)], (4) 
i 

where R; are the equilibrium positions of the charges, i. e. the underlying lattice sites. Using 
the Fourier representation of the delta function and solving (x) iteratively in terms of the 
displacement field, we obtain a decomposition of the density for the ith layer into a slowly 
and a rapidly spatially varying pieces: 

Pi (x)-n^-nV-u w (x)+ ]T ne iG,[x+uW(x)] , (5) 

G^O 

where G is a reciprocal lattice vector. Note that we have neglected terms that are products 
of the slowly and the rapidly varying terms. Physically, the first term represents density 
fluctuations for wavelengths greater than the lattice constant, and the second term represents 
the underlying lattice, modified by thermal fluctuations. Using the density decomposition 
(1), Hmt m &y be written as 

PH mt= f J^M^ e -i<i /> x /■ rf 2 xVq .(x-x') I nV . u A (x) _ y ne *G.[x + ^(x)]] 
J (2vr) 2 q J J y ^ Q J 

x ( nV-u B (x')- £ ne* G '^' +c+uB ^A , (6) 

V G'#0 / 

where c is the relative displacement vector between two lattices of the different plane and we 
have used the fact that ^ x l +d i = / A^s e* qx y e~ qd . Again neglecting the products of slowly 
and rapidly varying terms, which give vanishingly small contributions when integrating over 
all space, Hint separates into two pieces: a long-wavelength term 



and a short-wavelength term 



G^O G'^O 



(2vr) 2 g " C 



| rf 2 x | d 2 xV^ x - x ') e iG '[ x+uA Wl e i G'[ x ' +c+uS ( x ')]. (8) 



In order to obtain a tractable analytical treatment, we approximate this expression by 
splitting the sum over G' into two parts. The dominant part, with G' = — G is 

PHL = - E / 4^-2 ^-e-* d J d 2 X J dVe^ +G X X - X '> e ^.[u-(x)-u-(x')] j (Q) 



G^O 



where we have used e lG ' c = — 1. The second part (those terms with G' ^ — G) contains extra 
phase factors which tend to average to zero in the elastic limit. As a first approximation, 
we neglect such terms. Finally, Eq. (H) can be systematically expanded using a gradient 
expansion: 

mL = - E A G(d) f rf 2 x cos {G • [u A (x) - u B (x)] } + 0(d a ufd jU ^), (10) 

G^O J 

where Ag(gT) = 47r; g n2 ' e~ Gd . Putting Equations (P, (|7|), and ([K]) together, we obtain an 
effective Hamiltonian for the coupled planar Wigner crystals: 

P?U = PHo + J ^^^e-^ M/9 ^(q)«f(-q) 

- £ A G( d ) I rf ' x cos {G • [u A (x) - u B (x)] } . (11) 

G^O J 

The second term in Eq. ( O ) comes from the long- wavelength couplings while the third 
term reflects the periodicity of the underlying lattice structure. This particular structure 
in the effective Hamiltonian, as will be demonstrated below, leads to a total force which is 
comprised of two pieces - an exponentially decaying (short-ranged) force and a long-ranged 
power-law force: 

where Aq is the area of the plane. It is important to emphasize that both forces are present 
simultaneously, although each force dominates at a different spatial scale - the long-ranged 
force dominates at large separations while the short-ranged force at small separations. 



To calculate various expectation values in Eq. fll2|), it is convenient to transform the 
displacement fields into in-phase and out-of-phase displacement fields by u + (x) = u (x) + 
u B (x) and u~(x) = u A (x) — u B (x), respectively, so that the effective Hamiltonian ( [TT]) 
separates into two independent parts: 7i e = 7i + + H- with 

1 r d 2 q 



2 / (2^ n ^ (q) M ^ q ) M ^ _q )' 



(13) 



and 



1 .72 

m - = l\ f2^ n ^ (q) ^( q ) " £ / d 2 * cos[G ■ u-(x)], (14) 

where II^g(q) = 



G^O 



Ms^-(1 ± e -9 d ) + fiPap q 2 - Furthermore, at low temperature, 
where | u (x) | is small compared to the lattice constant a, the cosine term in Eq. ( $M can 
be expanded up to second order in |u~(x)| to obtain the "mass" terms. Within a harmonic 
approximation TC-, up to an additive constant, may be written as 
1 r d 2 Q „_ , , _, , , 1 f d\ 



d 2 q _ 1 



(27T) 



(15) 



Here, m 2 LT = AttIbu 2 Sg^o G e~ Gd = A^Ibu 2 A (d). Note that the mass terms vanish expo- 
nentially with d as also found in Ref. fl5[]. The fact that the transverse ttit and longitudinal 
"mass" m L are degenerate is related to the underlying triangular structure of the lattices 



| T5fl . These "masses" are associated with the finite energy required to uniformly shear the 
two Wigner crystals, and thus give rise to a gap in the dispersion relations of the out-of- 
phase modes. In the next two subsections, we derive expressions for the long-ranged and 
the short-ranged pressure as given in Eq. ( O) within the harmonic approximation. 



A. Long-ranged Pressure 



The long-ranged power-law force comes from the correlated long-wavelength density fluc- 
tuations (the plasmon modes). The shear modes do not contribute to this interaction since 
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d a PuB u a ' ( x ) = 0- Using Eqs. (^) and (|12[) , we obtain an expression for the long-ranged 
force: 



2nl, 



I 7^^(Wq)^B(-q)), 



(16) 



A J (2tt)2 

where 8pi(x) = — n V • u^(x) is the long-wavelength density fluctuation to the lowest order. 
Making use of the equipartition theorem to evaluate 



( 8 PA (q) 5p B (- q) } oc J Dii^q) 8p A (q) 5p B (- q) e~^ n + +n - 



-A, 



n 2 q 2 



Til B n 2 q (1 - e-Q d ) + ra\ 7il B n 2 q (1 + e~i d ) 



(17) 



and substituting this into Eq. (0), we have the result 

k B T 



d 3 



a(A d), 



(18) 



where 



, s C(3) x r 

a x - + _ r 



Ci(2v^) cos(2v^) + Si(2v/i) sin(2Vz) 



(19) 



C is the Riemann zeta function, and Ci(x) and Si (a;) are the cosine and sine integral functions 
TO , respectively. In the large distance limit, the second term in Eq. (|T9| ) is exponentially 



suppressed and can be neglected, yielding a = Therefore, for large d we have 



n LR (d) 



C(3) fc fl r 

8tt d 3 ' 



(20) 



This is the well-known result from the Debye-Hiickel approximation f7|. Note also that the 
amplitude = 0.048 is universal for this interaction, induced by the long wavelength 



fluctuations 17 . 



B. Short-ranged Pressure 



The short-ranged force which decays exponentially owes its existence to the "structural" 
correlations. It survives even at non-zero temperature, in contrast to the conclusion drawn 
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from a single 2D Wigner crystal, as discussed in the Introduction. However, we expect 
on physical grounds the short-ranged force to be weakened by thermal fluctuations. To 
compute its temperature dependence explicitly, we start with the expression for this force 
derived from Eqs. (0) and (|T2|): 



f3U SR (d) = -2vl B n 2 £ e ~— <l"~<°)l"> f G (d), (21) 

G^O 

where f G (d) = f ^S{q- G) e~ qd , 5(q - G) = J d2 re i(q-G)-r e -^[B+(r)-B-(r)] ) and 
S ± (r) = ( [u^r) — ir^O)] 2 ). Note that Eq. (|2~l"l) is exact, provided all the averages are 
evaluated exactly. For a system of coupled perfect Wigner crystals at zero temperature, 
/g(^) — e~ Gd - At finite temperature, but below the melting temperature T m , we note that 
B ± {r) varies very slowly in space, so that /g(^) can be approximated by its zero temperature 
value: fc,{d) ~ e~ Gd . Hence, we obtain 

pU SR (d) = -2nl B n 2 Y, e~ Gd (e^^W-^Wl) . (22 ) 

G^O 

The thermal average of the displacement fields in Eq. ( f22] ) resembles a "Debye- Waller" factor 
and indicate the degree to which the short-ranged force is depressed by thermal fluctuations 
from its zero temperature maximum value. Because of the cosine term present in Eq. (pT|) , 
this "Debye- Waller" factor is in general not zero, unlike the case of a single 2D Wigner 
crystal. However, if the system has melted into a Coulomb fluid, this cosine term, which 
comes from the lattice structure, would have to be modified. 

The required expectation value in Eq. fl2~2|) only involves 7i_. Within the harmonic 
approximation, the mean-square out-of-phase displacement field can be evaluated 



d 



AA (d)a 2 



+ In 

2nfi 



87rl B n 2 A (d)a 2 



God 
~2^ 



Xd 1 
nd fi 



(23) 



where Xr> — 1/ (2nl B n), a is the lattice constant, /x ~ 0.245 n 3//2 l B is the shear modulus of an 
isolated Wigner crystal in units of k B T, and in the last line, we have approximated Aq(c?) 
by the first nonzero reciprocal lattice vector contribution: Ao(d) ~ Go e~ G ° d . Note also that 
the logarithmic dependence on the "mass" (= Aim 2 l B Ao(d)) is a characteristic of 2D solids. 
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Inserting Eq. ( p3|) into Eq. (p^ ), we obtain an expression for the short-ranged pressure at 
finite temperatures 

/3U SR {d) ~ - 27r/ B n 2 e -( 1+ «/ 2 ) G ^. (24) 
Here, the parameter £ defined by 

characterizes the relative strengths of thermal fluctuations and the electrostatic energy of 
a Wigner crystal, i.e. £ ~ fcj ^J a - Thus, the sole effect of thermal fluctuations on the 
short-ranged force is to reduce its range: Go — > Go ( 1 + f 

C. Discussion of Results 

In summary, we have shown that the total pressure can be decomposed into a long-ranged 
II and a short-ranged pressure Hsr- Each force is computed at low temperatures, where 
the harmonic approximation is expected to be valid. The result for the total force is 

(311(d) ^-2nl B n 2 e-^ + ^ God - 2<M f ( 26 ) 

(Jj 

C(3) 



where £ = ^ ^4 + - and a(A <i) = for large d. In Fig. §, we have plotted Hsr and 



Ad i 1 



nd fj, 



Ulr for two values of the coupling constant, r = — 150 and 50. Not surprisingly, they 
show that Hsr dominates for small d, and for large d. However, it is interesting to note 
that even for high values of T, Ulr dominates as soon as d ~ a. 



According to Eq. (|24|), the magnitude of Hsr tends to decrease exponentially with 
temperature, as illustrated in Fig. [3| This strong decrease with increasing temperature 
is consistent with the Brownian dynamics simulations of Gr0nbech- Jensen et al. 0. The 
shortening of its range may be attributed to the generic nature of strong fluctuations in 2D 
systems, and can also be understood by the following scaling argument. Referring back to 

in Eq. flUD, one can show that the anomalous dimension of the operator cos[Go • u~(x)] 
is [ Length and correspondingly the dimension of Ag (g?) is [ Length ]^~ 2 . Since Aq (c?) 
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is the only relevant length scale in 7Y_, we must have (V Go ' u (°n ~ A G 2 Q c JT^|. Therefore, 
the short-ranged pressure scales like 

U SR (d) ~ — Ag (d) x (e iG °' u -(°)) ~ - A Go x A G f* ~ - e ' God W. (27) 



J l-l Gq ^ ^ Gq 

In the low temperature limit < 1), we see that the range of Hsr is G ( 1 + § ) as * n Eq.- 



(24j). This scaling argument also suggests that at higher temperatures thermal fluctuations 
may have interesting nonperturbative effects. At zero temperature £ = 0, so in Eq. 
(f2~4|) reproduces the known result of exponentially decaying attractive force |||§. 

The long-ranged pressure for large d in Eq. (|20| ) agrees exactly, including the prefactor, 
with the Debye-Hiickel approximation. This is hardly surprising since the existence of long- 
wavelength plasmons (average density fluctuations) is independent of local structure, and 
they are present for solids and fluids alike. Thus, the asymptotic long-ranged power-law 
force must manifest itself even after QLRO is lost via a 2D melting transition driven by 



dislocations ||19|| . Therefore, our formulation captures the essential physics of the attraction 
not only arising from the ground state "structural" correlations, but also from the high 
temperature charge-fluctuations. 



III. QUANTUM CONTRIBUTIONS TO THE LONG-RANGED ATTRACTION 

According to the classical calculations above, correlation effects give rise to a "struc- 
tural" short-ranged and a long-ranged attractive force. Recall that the long-ranged force 
vanishes as T — > 0, and that the short-ranged force is strongest at zero temperature but 
vanishes exponentially with distance. This observation suggests that for sufficiently large 
separations correlated attractions at finite temperatures are stronger than those arising from 
the perfectly correlated zero temperature ground state. However, we have pointed out in 
Ref. H that zero-point fluctuations of the plasmons lead to an attractive long-ranged in- 
teraction, which exhibits an unusual fractional-power-law decay (~ d™ 7 / 2 ), in contrast to 
the zero-temperature van der Waals interaction (~ d -4 ). Hence, in the T — ► limit, this 
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"zero-point attraction" dominates the short-ranged "structural" force at large separations. 
Furthermore, we expect that quantum fluctuations persist at finite temperature, and in this 
section, we compute their temperature dependence. 

Within the harmonic approximation to the effective Hamiltonian, the dispersion relations 
for the plasmons can be readily obtained ||: 

u?M = — A (d) + — g ( 1 - e~ qd ); (28) 
ul( q ) = 2 -^ q (l + e^ d ), (29) 

where m is the mass of the charges and A (d) ~ e~ Gd is proportional to the energy gap (the 
"mass" term) for the out-of-phase mode. The plasmon modes are related to the correlated 
charge- density fluctuations in the two layers. At any finite temperature, the free energy of 
the low-lying plasmon excitations is given by 

mi A, = \ E / (g, «.(q) H- k B T £ / *« ln (1 _ , (30) 

where A Q is the area of the plane. Since the energy gap A is exponentially damped for 
large distances, its contribution to the free energy may be neglected in the large distance 
limit, where the long-ranged force is expected to be dominant. 

The first term in Eq. ( |3T)D arising from the zero-point fluctuations has been computed in 
Ref. @ and gives the d~ 7 ^ 2 power-law mentioned above. An additional contribution to the 
pressure at finite temperature arises from the second term in Eq. (|30|), 

HA f°° . K/n i 1 



^L R {d) = --^ m HdxX^l y= 

4ndV 2 Jo [exp[riJx(l-e-*)]-l Vl 



1 



exp[ 



[V\/x(l + e~ x )} - 1 Vl + e~ x 



(31) 



where A = and V = PhA/Vd. We can evaluate this expression in two limits: 

In the low-temperature limit rj ^> 1, Eq. ([!!]) can be systematically expanded in powers 
of ?7 _1 . The lowest order term is given by Hui(d) = —®j^j2, where = a-Bjx^, a B = 
eh 2 j (me 2 ) is the effective Bohr radius, a = j- f °° dx -§^j = C(3)/ (2vr), and ( is the Riemann 
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zeta function. We observe that the low temperature condition 77 > 1 is equivalent to the 
short distance limit d < Xl- 

In the high temperature limit 77 <C 1 or the large distance limit d > Xl, we expand the ex- 
ponential in the denominator of Eq. fl3T|) to obtain TlLR(d) = —a -J^, where a = C(3)/(87r). 
This result agrees with the classical calculation in Sec. [11 A| as it should. Therefore, we have 
the following regimes for correlated attraction from plasmon fluctuations at finite tempera- 
ture 

-k B T/d 3 , for Xl < d, 
IL LR (d)~{ (32) 

[ -k B T/(X L d 2 ), for X L > d. 

We note that A^, in contrast to A^, increases with decreasing temperature, indicating, as one 
might expect, that quantum fluctuations are important at low temperatures. Furthermore, 
since IL^d) — > as T — > 0, the attractive interaction as T — > is governed by zero- 
point fluctuations as emphasized above. In the strong Coulomb coupling limit Ib/Xd ~ 
100, we get Xl ~ 3Afor e ~ 100 and a B ~ 1/20A. Finally, it should be emphasized 
that quantum contributions to the long-ranged attraction are unlikely to be relevant for 
macroions. Our motivation here stems from the desire to understand the charge-fluctuation- 
induced attraction between coupled layers in a complete picture. However, our results 
may be relevant for electrons in bilayer semiconductor systems. Indeed, there are recent 
theoretical efforts devoted to this subject [ED . 



IV. DISCUSSION AND CONCLUSION 

In this paper, we have studied analytically the electrostatic attraction between two planar 
Wigner crystals in the strong Coulomb coupling limit. We show that the total attractive 
pressure can be separated into a long-ranged and short-ranged component. The long-ranged 
pressure arises from correlated fluctuations and the short-ranged pressure from the ground 
state "structural" correlations. We also compute the very low temperature behavior of the 
fluctuation-induced attraction, where long-wavelength plasmon excitation must be described 
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by Bose-Einstein statistics. The results are summarized in Fig. |], showing different regimes 
for the charge-fluctuation-induced long-ranged attraction, including the high temperature 
results in Ref. and the characteristic decay length 1$r for the short-ranged force. For small 
d, the short-ranged force is always dominant, but the decay length shrinks with increasing 
temperature. The crossover from the short-ranged to long-ranged dominant regimes occurs 
about d ~ a. Thus, for large d 3> a only the long-ranged force is operative, which crosses 
over from d~ 7 ^ 2 at zero temperature to the finite temperature distance dependence of d~ 2 if 
d < \l and d~ 3 if d > Xl- This provides a unified description to the electrostatic attraction 
between two coupled Wigner crystals. 

In addition, our formulation may offer further insights into the nature of the counterion- 



mediated attraction at short distances. As discussed in Sec. [II B| , the reason that the short- 
ranged force in Eq. (p^ ) does not vanish is because of the cosine term in Ti-, which represents 
the underlying lattice structures, and our results indicate that the strength of the short- 
ranged force decreases exponentially with temperature. However, at higher temperatures 
the expression for Hsr in Eq. (^) is no longer valid, since the harmonic approximation 
breaks down. Indeed, the scaling argument leading to Eq. ( P7| ) suggests that if the full 
cosine term is retained, Hsr may exhibit nonperturbative behaviors as £ — > 2~. 

To discuss qualitatively what happens at higher temperatures, we assume that Aa(d) is 
sufficiently small and the system of interacting Wigner crystals is below its melting temper- 
ature T m . Then, the charges between the two layers may unlock via a Kosterlitz-Thouless 
(KT) type of transition, determined by the relevancy of the cosine term in 7Y_, at £ = 2 
pl| . (An order of magnitude estimate for the coupling constant is T ~ 13.) In the locked 
phase, ^ <C 2, the periodic symmetry in H- is spontaneously broken, and the resulting 
state is well captured by the harmonic approximation. On the other hand, when £ > 2 the 
fluctuations are so large that the ground state becomes nondegenerate (gapless), i.e. the 
layers are decoupled. To compute Hsr in the unlocked phase, 7if nt given in Eq. ( |i0|) can 
be treated as a perturbation in evaluating the "Debye- Waller" factor in Eq. (E3). To the 



lowest order, we obtain 
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n -<^-fHf^ )e " (33) 



We first note that this expression diverges as £ — > 2 + , indicating the breakdown of the 
perturbation theory as the temperature is lowered. Furthermore, in contrast to Eq. (p4"|), 
the range of Hsr remains constant and the amplitude acquires a temperature dependence 
of ~ 1/T (for large £ ^> 2), reminiscent of a high temperature expansion. 

However, the above picture may be modified if the charges have melted into a Coulomb 
fluid via a dislocation-mediated melting transition |L9[ before £ —>■ 2~. If this is the case, 
further analysis is necessary to obtain a more complete picture of the high temperature 
phase. Although the spatial correlations in a system of coupled 2D Coulomb fluids are 
expected to be somewhat different from 2D Wigner crystals, the solid phase results above 
suggest a qualitative lower limit of T ~ 13 at which Hsr crosses over from low temperature 
to high temperature behavior. It may be of interest to note that in Ref. |J, an estimate for 
the upper limit of T at which the Poisson-Boltzmann equation breaks down is of the order 
of r ~ 3. To describe the melting of coupled 2D Wigner crystals, excitations of dislocations 
must be introduced into the effective Hamiltonian Eq. ( |TTD similar to what is done in Ref. 



1 22]] . These considerations may help to establish an analytical theory of the attraction arising 
from counterion correlations, not captured by the Poisson-Boltzmann theory. The present 
formulation is a first step in that direction. 
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FIGURES 





FIG. 1. A schematic picture of two staggered Wigner crystals formed by the "condensed : 
counterions at very low temperatures. 
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FIG. 3. A plot of Hsr as a function of T 1 at d = 1.2 a, according to Eq. (|24|), which shows 
that TIsr exponentially decreases with increasing temperature. 




FIG. 4. A schematic phase diagram summarizing different charge-fiuctuation-induced attrac- 
tion regimes. The characteristic decay length Isr of the short-ranged force is also shown. 
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